home *** CD-ROM | disk | FTP | other *** search
/ Ham Radio 2000 / Ham Radio 2000.iso / ham2000 / misc / dspice0s / noise.c < prev    next >
C/C++ Source or Header  |  1992-11-21  |  48KB  |  1,446 lines
  1. /* noise.f -- translated by f2c (version of 3 February 1990  3:36:42).
  2.    You must link the resulting object file with the libraries:
  3.     -lF77 -lI77 -lm -lc   (in that order)
  4. */
  5.  
  6. #include "f2c.h"
  7.  
  8. /* Common Block Declarations */
  9.  
  10. struct {
  11.     integer ielmnt, isbckt, nsbckt, iunsat, nunsat, itemps, numtem, isens, 
  12.         nsens, ifour, nfour, ifield, icode, idelim, icolum, insize, 
  13.         junode, lsbkpt, numbkp, iorder, jmnode, iur, iuc, ilc, ilr, 
  14.         numoff, isr, nmoffc, iseq, iseq1, neqn, nodevs, ndiag, iswap, 
  15.         iequa, macins, lvnim1, lx0, lvn, lynl, lyu, lyl, lx1, lx2, lx3, 
  16.         lx4, lx5, lx6, lx7, ld0, ld1, ltd, imynl, imvn, lcvn, nsnod, 
  17.         nsmat, nsval, icnod, icmat, icval, loutpt, lpol, lzer, irswpf, 
  18.         irswpr, icswpf, icswpr, irpt, jcpt, irowno, jcolno, nttbr, nttar, 
  19.         lvntmp;
  20. } tabinf_;
  21.  
  22. #define tabinf_1 tabinf_
  23.  
  24. struct {
  25.     integer locate[50], jelcnt[50], nunods, ncnods, numnod, nstop, nut, nlt, 
  26.         nxtrm, ndist, ntlin, ibr, numvs, numalt, numcyc;
  27. } cirdat_;
  28.  
  29. #define cirdat_1 cirdat_
  30.  
  31. struct {
  32.     doublereal omega, time, delta, delold[7], ag[7], vt, xni, egfet, xmu, 
  33.         sfactr;
  34.     integer mode, modedc, icalc, initf, method, iord, maxord, noncon, iterno, 
  35.         itemno, nosolv, modac, ipiv, ivmflg, ipostp, iscrch, iofile;
  36. } status_;
  37.  
  38. #define status_1 status_
  39.  
  40. struct {
  41.     doublereal atime, aprog[3], adate, atitle[10], defl, defw, defad, defas, 
  42.         rstats[50];
  43.     integer iwidth, lwidth, nopage;
  44. } miscel_;
  45.  
  46. #define miscel_1 miscel_
  47.  
  48. struct {
  49.     doublereal twopi, xlog2, xlog10, root2, rad, boltz, charge, ctok, gmin, 
  50.         reltol, abstol, vntol, trtol, chgtol, eps0, epssil, epsox, pivtol,
  51.          pivrel;
  52. } knstnt_;
  53.  
  54. #define knstnt_1 knstnt_
  55.  
  56. struct {
  57.     doublereal fstart, fstop, fincr, skw2, refprl, spw2;
  58.     integer jacflg, idfreq, inoise, nosprt, nosout, nosin, idist, idprt;
  59. } ac_;
  60.  
  61. #define ac_1 ac_
  62.  
  63. struct {
  64.     doublereal value[200000];
  65. } blank_;
  66.  
  67. #define blank_1 blank_
  68.  
  69. /* Table of constant values */
  70.  
  71. static integer c__1 = 1;
  72. static integer c__2 = 2;
  73. static integer c__0 = 0;
  74. static integer c__8 = 8;
  75. static integer c__5 = 5;
  76.  
  77. /*<       subroutine noise(loco) >*/
  78. /* Subroutine */ int noise_(loco)
  79. integer *loco;
  80. {
  81.     /* Initialized data */
  82.  
  83.     static struct {
  84.     char e_1[8];
  85.     doublereal e_2;
  86.     } equiv_106 = { {'r', 'c', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  87.  
  88. #define alsrc (*(doublereal *)&equiv_106)
  89.  
  90.     static struct {
  91.     char e_1[8];
  92.     doublereal e_2;
  93.     } equiv_107 = { {'r', 'e', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  94.  
  95. #define alsre (*(doublereal *)&equiv_107)
  96.  
  97.     static struct {
  98.     char e_1[8];
  99.     doublereal e_2;
  100.     } equiv_108 = { {'r', 's', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  101.  
  102. #define alsrs (*(doublereal *)&equiv_108)
  103.  
  104.     static struct {
  105.     char e_1[8];
  106.     doublereal e_2;
  107.     } equiv_109 = { {'r', 'd', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  108.  
  109. #define alsrd (*(doublereal *)&equiv_109)
  110.  
  111.     static struct {
  112.     char e_1[8];
  113.     doublereal e_2;
  114.     } equiv_110 = { {'i', 'b', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  115.  
  116. #define alsib (*(doublereal *)&equiv_110)
  117.  
  118.     static struct {
  119.     char e_1[8];
  120.     doublereal e_2;
  121.     } equiv_111 = { {'i', 'c', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  122.  
  123. #define alsic (*(doublereal *)&equiv_111)
  124.  
  125.     static struct {
  126.     char e_1[8];
  127.     doublereal e_2;
  128.     } equiv_112 = { {'i', 'd', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  129.  
  130. #define alsid (*(doublereal *)&equiv_112)
  131.  
  132.     static struct {
  133.     char e_1[8];
  134.     doublereal e_2;
  135.     } equiv_113 = { {'f', 'n', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  136.  
  137. #define alsfn (*(doublereal *)&equiv_113)
  138.  
  139.     static struct {
  140.     char e_1[8];
  141.     doublereal e_2;
  142.     } equiv_114 = { {'t', 'o', 't', 'a', 'l', ' ', ' ', ' '}, 0. };
  143.  
  144. #define alstot (*(doublereal *)&equiv_114)
  145.  
  146.     static struct {
  147.     char e_1[8];
  148.     doublereal e_2;
  149.     } equiv_115 = { {'/', ' ', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  150.  
  151. #define aslash (*(doublereal *)&equiv_115)
  152.  
  153.     static struct {
  154.     char e_1[8];
  155.     doublereal e_2;
  156.     } equiv_116 = { {' ', ' ', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  157.  
  158. #define ablnk (*(doublereal *)&equiv_116)
  159.  
  160.     static struct {
  161.     char e_1[24];
  162.     doublereal e_2;
  163.     } equiv_117 = { {'(', '/', '/', '/', '/', ',', '1', '1', 'x', ',', ' '
  164.         , ' ', '(', '2', 'x', ',', 'a', '8', ')', ')', ' ', ' ', ' ', 
  165.         ' '}, 0. };
  166.  
  167. #define afmt1 ((doublereal *)&equiv_117)
  168.  
  169.     static struct {
  170.     char e_1[24];
  171.     doublereal e_2;
  172.     } equiv_118 = { {'(', '1', 'h', '0', ',', 'a', '8', ',', '1', 'p', ' '
  173.         , ' ', 'd', '1', '0', '.', '3', ')', ' ', ' ', ' ', ' ', ' ', 
  174.         ' '}, 0. };
  175.  
  176. #define afmt2 ((doublereal *)&equiv_118)
  177.  
  178.     static struct {
  179.     char e_1[32];
  180.     doublereal e_2;
  181.     } equiv_119 = { {'n', 'o', 'i', 's', 'e', ' ', 'a', 'n', 'a', 'l', 'y'
  182.         , 's', 'i', 's', ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', 
  183.         ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  184.  
  185. #define titln ((doublereal *)&equiv_119)
  186.  
  187.     static struct {
  188.     char e_1[8];
  189.     doublereal e_2;
  190.     } equiv_120 = { {'r', 'b', ' ', ' ', ' ', ' ', ' ', ' '}, 0. };
  191.  
  192. #define alsrb (*(doublereal *)&equiv_120)
  193.  
  194.  
  195.     /* Format strings */
  196.     static char fmt_16[] = "(\0020    frequency = \002,1pd10.3,\002 hz\002/)";
  197.  
  198.     static char fmt_91[] = "(//\0020**** resistor squared noise voltages (sq\
  199.  v/hz)\002)";
  200.     static char fmt_201[] = "(//\0020**** diode squared noise voltages (sq v\
  201. /hz)\002)";
  202.     static char fmt_301[] = "(//\0020**** transistor squared noise voltages \
  203. (sq v/hz)\002)";
  204.     static char fmt_401[] = "(//\0020**** jfet squared noise voltages (sq v/\
  205. hz)\002)";
  206.     static char fmt_501[] = "(//\0020**** mosfet squared noise voltages (sq \
  207. v/hz)\002)";
  208.     static char fmt_611[] = "(////,\0020**** total output noise voltage\002,\
  209. 9x,\002= \002,1pd10.3,\002 sq v/hz\002/,\0020\002,40x,\002= \002,d10.3,\002 \
  210. v/rt hz\002/,\0020     transfer function value:\002,/,\0020\002,7x,4a8,a1\
  211. ,\002= \002,d10.3,/,\0020     equivalent input noise at \002,a8,\002 = \002,\
  212. d10.3,\002 /rt hz\002)";
  213.  
  214.     /* System generated locals */
  215.     integer i_1, i_2;
  216.     doublereal d_1, d_2;
  217.     complex q_1;
  218.     static complex equiv_3[1];
  219.  
  220.     /* Builtin functions */
  221.     integer s_wsfe(), do_fio(), e_wsfe();
  222.     double sqrt(), log(), exp();
  223.  
  224.     /* Local variables */
  225.     static doublereal area, anam[12];
  226. #define cval (equiv_3)
  227.     static integer locm;
  228.     static doublereal freq;
  229.     static integer locv;
  230.     extern /* Subroutine */ int move_();
  231.     static integer ipos, loct, kntr;
  232.     static doublereal gdpr, gspr, vnin, twoq, vout, anam1;
  233.     static integer node1, node2, node3, node4, node5, node6;
  234. #define c (equiv_3)
  235.     extern /* Subroutine */ int zero8_();
  236.     static integer i;
  237. #define v ((real *)equiv_3)
  238. #define vimag ((real *)equiv_3 + 1)
  239. #define vreal ((real *)equiv_3)
  240.     static doublereal vntot[12];
  241.     static integer nprnt, nkntr;
  242.     extern /* Subroutine */ int title_();
  243.     static doublereal vnrms, vtemp;
  244.     static integer itype;
  245.     static doublereal vnout, gm;
  246.     extern /* Subroutine */ int acasol_();
  247.     static doublereal xl;
  248. #define nodplc ((integer *)&blank_1)
  249. #define cvalue ((complex *)&blank_1)
  250.     static doublereal string[5];
  251.     extern /* Subroutine */ int alfnum_();
  252.     static doublereal fourkt;
  253.     static integer noposo, nonego, kntlim, ititle;
  254.     static doublereal fna;
  255.     static integer loc;
  256.     static doublereal fnk, arg, xw;
  257.     static integer ioutyp;
  258.     extern /* Subroutine */ int outnam_();
  259.     static doublereal cox, vno1[12], vno2[12], vno3[12], vno4[12], vno5[12], 
  260.         vno6[12];
  261.  
  262.     /* Fortran I/O blocks */
  263.     static cilist io__33 = { 0, 0, 0, fmt_16, 0 };
  264.     static cilist io__43 = { 0, 0, 0, fmt_91, 0 };
  265.     static cilist io__44 = { 0, 0, 0, (char *)afmt1, 0 };
  266.     static cilist io__46 = { 0, 0, 0, (char *)afmt2, 0 };
  267.     static cilist io__58 = { 0, 0, 0, fmt_201, 0 };
  268.     static cilist io__59 = { 0, 0, 0, (char *)afmt1, 0 };
  269.     static cilist io__60 = { 0, 0, 0, (char *)afmt2, 0 };
  270.     static cilist io__61 = { 0, 0, 0, (char *)afmt2, 0 };
  271.     static cilist io__62 = { 0, 0, 0, (char *)afmt2, 0 };
  272.     static cilist io__63 = { 0, 0, 0, (char *)afmt2, 0 };
  273.     static cilist io__70 = { 0, 0, 0, fmt_301, 0 };
  274.     static cilist io__71 = { 0, 0, 0, (char *)afmt1, 0 };
  275.     static cilist io__72 = { 0, 0, 0, (char *)afmt2, 0 };
  276.     static cilist io__73 = { 0, 0, 0, (char *)afmt2, 0 };
  277.     static cilist io__74 = { 0, 0, 0, (char *)afmt2, 0 };
  278.     static cilist io__75 = { 0, 0, 0, (char *)afmt2, 0 };
  279.     static cilist io__76 = { 0, 0, 0, (char *)afmt2, 0 };
  280.     static cilist io__77 = { 0, 0, 0, (char *)afmt2, 0 };
  281.     static cilist io__78 = { 0, 0, 0, (char *)afmt2, 0 };
  282.     static cilist io__79 = { 0, 0, 0, fmt_401, 0 };
  283.     static cilist io__80 = { 0, 0, 0, (char *)afmt1, 0 };
  284.     static cilist io__81 = { 0, 0, 0, (char *)afmt2, 0 };
  285.     static cilist io__82 = { 0, 0, 0, (char *)afmt2, 0 };
  286.     static cilist io__83 = { 0, 0, 0, (char *)afmt2, 0 };
  287.     static cilist io__84 = { 0, 0, 0, (char *)afmt2, 0 };
  288.     static cilist io__85 = { 0, 0, 0, (char *)afmt2, 0 };
  289.     static cilist io__93 = { 0, 0, 0, fmt_501, 0 };
  290.     static cilist io__94 = { 0, 0, 0, (char *)afmt1, 0 };
  291.     static cilist io__95 = { 0, 0, 0, (char *)afmt2, 0 };
  292.     static cilist io__96 = { 0, 0, 0, (char *)afmt2, 0 };
  293.     static cilist io__97 = { 0, 0, 0, (char *)afmt2, 0 };
  294.     static cilist io__98 = { 0, 0, 0, (char *)afmt2, 0 };
  295.     static cilist io__99 = { 0, 0, 0, (char *)afmt2, 0 };
  296.     static cilist io__105 = { 0, 0, 0, fmt_611, 0 };
  297.  
  298.  
  299. /*<       implicit double precision (a-h,o-z) >*/
  300.  
  301. /*     this routine computes the noise due to various circuit elements. */
  302.  
  303.  
  304. /* spice version 2g.6  sccsid=tabinf 3/15/83 */
  305. /*<       common /tabinf/ ielmnt,isbckt,nsbckt,iunsat,nunsat,itemps,numtem, >*/
  306. /*<      1   isens,nsens,ifour,nfour,ifield,icode,idelim,icolum,insize, >*/
  307. /*<      2   junode,lsbkpt,numbkp,iorder,jmnode,iur,iuc,ilc,ilr,numoff,isr, >*/
  308. /*<      3   nmoffc,iseq,iseq1,neqn,nodevs,ndiag,iswap,iequa,macins,lvnim1, >*/
  309. /*<      4   lx0,lvn,lynl,lyu,lyl,lx1,lx2,lx3,lx4,lx5,lx6,lx7,ld0,ld1,ltd, >*/
  310. /*<      5   imynl,imvn,lcvn,nsnod,nsmat,nsval,icnod,icmat,icval, >*/
  311. /*<      6   loutpt,lpol,lzer,irswpf,irswpr,icswpf,icswpr,irpt,jcpt, >*/
  312. /*<      7   irowno,jcolno,nttbr,nttar,lvntmp >*/
  313. /* spice version 2g.6  sccsid=cirdat 3/15/83 */
  314. /*<       common /cirdat/ locate(50),jelcnt(50),nunods,ncnods,numnod,nstop, >*/
  315. /*<      1   nut,nlt,nxtrm,ndist,ntlin,ibr,numvs,numalt,numcyc >*/
  316. /* spice version 2g.6  sccsid=status 3/15/83 */
  317. /*<       common /status/ omega,time,delta,delold(7),ag(7),vt,xni,egfet, >*/
  318. /*<      1   xmu,sfactr,mode,modedc,icalc,initf,method,iord,maxord,noncon, >*/
  319. /*<      2   iterno,itemno,nosolv,modac,ipiv,ivmflg,ipostp,iscrch,iofile >*/
  320. /* spice version 2g.6  sccsid=miscel 3/15/83 */
  321. /*<       common /miscel/ atime,aprog(3),adate,atitle(10),defl,defw,defad, >*/
  322. /*<      1  defas,rstats(50),iwidth,lwidth,nopage >*/
  323. /* spice version 2g.6  sccsid=knstnt 3/15/83 */
  324. /*<       common /knstnt/ twopi,xlog2,xlog10,root2,rad,boltz,charge,ctok, >*/
  325. /*<      1   gmin,reltol,abstol,vntol,trtol,chgtol,eps0,epssil,epsox, >*/
  326. /*<      2   pivtol,pivrel >*/
  327. /* spice version 2g.6  sccsid=ac 3/15/83 */
  328. /*<       common /ac/ fstart,fstop,fincr,skw2,refprl,spw2,jacflg,idfreq, >*/
  329. /*<      1   inoise,nosprt,nosout,nosin,idist,idprt >*/
  330. /* spice version 2g.6  sccsid=blank 3/15/83 */
  331. /*<       common /blank/ value(200000) >*/
  332. /*<       integer nodplc(64) >*/
  333. /*<       complex cvalue(32) >*/
  334. /*<       equivalence (value(1),nodplc(1),cvalue(1)) >*/
  335.  
  336.  
  337. /*<       dimension vno1(12),vno2(12),vno3(12),vno4(12),vno5(12),vno6(12) >*/
  338. /*<       dimension vntot(12),anam(12),string(5) >*/
  339. /*<       real v,vreal,vimag >*/
  340. /*<       dimension titln(4),v(2) >*/
  341. /*<       dimension afmt1(3),afmt2(3) >*/
  342. /*<       complex cval,c(1) >*/
  343. /*<       equivalence (c(1),v(1),cval) >*/
  344. /*<       equivalence (v(1),vreal),(v(2),vimag) >*/
  345. /*<       data titln / 8hnoise an, 8halysis  , 8h        , 8h         / >*/
  346. /*<       data alsrb,alsrc,alsre,alsrs,alsrd / 2hrb,2hrc,2hre,2hrs,2hrd / >*/
  347. /*<       data alsib,alsic,alsid,alsfn / 2hib,2hic,2hid,2hfn / >*/
  348. /*<       data alstot / 5htotal / >*/
  349. /*<       data aslash,ablnk / 1h/, 1h  / >*/
  350. /*<       data afmt1 /8h(////,11,8hx,  (2x,,8ha8))    / >*/
  351. /*<       data afmt2 /8h(1h0,a8,,8h1p  d10.,8h3)      / >*/
  352.  
  353.  
  354. /* .. fix-up formats */
  355. /*<       kntr=12 >*/
  356.     kntr = 12;
  357. /*<       if(lwidth.le.80) kntr=7 >*/
  358.     if (miscel_1.lwidth <= 80) {
  359.     kntr = 7;
  360.     }
  361. /*<       ipos=11 >*/
  362.     ipos = 11;
  363. /*<       call move(afmt1,ipos,ablnk,1,2) >*/
  364.     move_(afmt1, &ipos, &ablnk, &c__1, &c__2);
  365. /*<       call alfnum(kntr,afmt1,ipos) >*/
  366.     alfnum_(&kntr, afmt1, &ipos);
  367. /*<       ipos=11 >*/
  368.     ipos = 11;
  369. /*<       call move(afmt2,ipos,ablnk,1,2) >*/
  370.     move_(afmt2, &ipos, &ablnk, &c__1, &c__2);
  371. /*<       call alfnum(kntr,afmt2,ipos) >*/
  372.     alfnum_(&kntr, afmt2, &ipos);
  373. /*<       nprnt=0 >*/
  374.     nprnt = 0;
  375. /*<       freq=omega/twopi >*/
  376.     freq = status_1.omega / knstnt_1.twopi;
  377. /*<       if (icalc.ge.2) go to 10 >*/
  378.     if (status_1.icalc >= 2) {
  379.     goto L10;
  380.     }
  381. /*<       fourkt=4.0d0*charge*vt >*/
  382.     fourkt = knstnt_1.charge * 4. * status_1.vt;
  383. /*<       twoq=2.0d0*charge >*/
  384.     twoq = knstnt_1.charge * 2.;
  385. /*<       noposo=nodplc(nosout+2) >*/
  386.     noposo = nodplc[ac_1.nosout + 1];
  387. /*<       nonego=nodplc(nosout+3) >*/
  388.     nonego = nodplc[ac_1.nosout + 2];
  389. /*<       kntlim=lwidth/11 >*/
  390.     kntlim = miscel_1.lwidth / 11;
  391. /*<       nkntr=1 >*/
  392.     nkntr = 1;
  393. /*<    10 if (nosprt.eq.0) go to 30 >*/
  394. L10:
  395.     if (ac_1.nosprt == 0) {
  396.     goto L30;
  397.     }
  398. /*<       if (nkntr.gt.icalc) go to 30 >*/
  399.     if (nkntr > status_1.icalc) {
  400.     goto L30;
  401.     }
  402. /*<       nprnt=1 >*/
  403.     nprnt = 1;
  404. /*<       nkntr=nkntr+nosprt >*/
  405.     nkntr += ac_1.nosprt;
  406. /*<       call title(0,lwidth,1,titln) >*/
  407.     title_(&c__0, &miscel_1.lwidth, &c__1, titln);
  408. /*<       write (iofile,16) freq >*/
  409.     io__33.ciunit = status_1.iofile;
  410.     s_wsfe(&io__33);
  411.     do_fio(&c__1, (char *)&freq, (ftnlen)sizeof(doublereal));
  412.     e_wsfe();
  413. /*<    16 format('0    frequency = ',1pd10.3,' hz'/) >*/
  414.  
  415. /*  obtain adjoint circuit solution */
  416.  
  417. /*<    30 vnrms=0.0d0 >*/
  418. L30:
  419.     vnrms = 0.;
  420. /*<       cval=cvalue(lcvn+noposo)-cvalue(lcvn+nonego) >*/
  421.     i_1 = tabinf_1.lcvn + noposo - 1;
  422.     i_2 = tabinf_1.lcvn + nonego - 1;
  423.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  424.         .i;
  425.     cval->r = q_1.r, cval->i = q_1.i;
  426. /*<       vout=dsqrt(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  427.     vout = sqrt((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * *
  428.         vimag));
  429. /*<       vout=dmax1(vout,1.0d-20) >*/
  430.     vout = max(vout,1e-20);
  431. /*<       call zero8(value(lvn+1),nstop) >*/
  432.     zero8_(&blank_1.value[tabinf_1.lvn], &cirdat_1.nstop);
  433. /*<       call zero8(value(imvn+1),nstop) >*/
  434.     zero8_(&blank_1.value[tabinf_1.imvn], &cirdat_1.nstop);
  435. /*<       value(lvn+noposo)=-1.0d0 >*/
  436.     blank_1.value[tabinf_1.lvn + noposo - 1] = -1.;
  437. /*<       value(lvn+nonego)=+1.0d0 >*/
  438.     blank_1.value[tabinf_1.lvn + nonego - 1] = 1.;
  439. /*<       call acasol >*/
  440.     acasol_();
  441.  
  442. /*  resistors */
  443.  
  444. /*<       if (jelcnt(1).eq.0) go to 200 >*/
  445.     if (cirdat_1.jelcnt[0] == 0) {
  446.     goto L200;
  447.     }
  448. /*<       ititle=0 >*/
  449.     ititle = 0;
  450. /*<    91 format(//'0**** resistor squared noise voltages (sq v/hz)') >*/
  451. /* L91: */
  452. /*<   100 loc=locate(1) >*/
  453. /* L100: */
  454.     loc = cirdat_1.locate[0];
  455. /*<       kntr=0 >*/
  456.     kntr = 0;
  457. /*<   110 if ((loc.eq.0).or.(nodplc(loc+8).ne.0)) go to 130 >*/
  458. L110:
  459.     if (loc == 0 || nodplc[loc + 7] != 0) {
  460.     goto L130;
  461.     }
  462. /*<       kntr=kntr+1 >*/
  463.     ++kntr;
  464. /*<       locv=nodplc(loc+1) >*/
  465.     locv = nodplc[loc];
  466. /*<       anam(kntr)=value(locv) >*/
  467.     anam[kntr - 1] = blank_1.value[locv - 1];
  468. /*<       node1=nodplc(loc+2) >*/
  469.     node1 = nodplc[loc + 1];
  470. /*<       node2=nodplc(loc+3) >*/
  471.     node2 = nodplc[loc + 2];
  472. /*<       cval=cvalue(lcvn+node1)-cvalue(lcvn+node2) >*/
  473.     i_1 = tabinf_1.lcvn + node1 - 1;
  474.     i_2 = tabinf_1.lcvn + node2 - 1;
  475.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  476.         .i;
  477.     cval->r = q_1.r, cval->i = q_1.i;
  478. /*<       vntot(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  479. /*<      1   *fourkt*value(locv+1) >*/
  480.     vntot[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag *
  481.          *vimag)) * fourkt * blank_1.value[locv];
  482. /*<       vnrms=vnrms+vntot(kntr) >*/
  483.     vnrms += vntot[kntr - 1];
  484. /*<       if (kntr.ge.kntlim) go to 140 >*/
  485.     if (kntr >= kntlim) {
  486.     goto L140;
  487.     }
  488. /*<   120 loc=nodplc(loc) >*/
  489. L120:
  490.     loc = nodplc[loc - 1];
  491. /*<       go to 110 >*/
  492.     goto L110;
  493. /*<   130 if (kntr.eq.0) go to 200 >*/
  494. L130:
  495.     if (kntr == 0) {
  496.     goto L200;
  497.     }
  498. /*<   140 if (nprnt.eq.0) go to 160 >*/
  499. L140:
  500.     if (nprnt == 0) {
  501.     goto L160;
  502.     }
  503. /*<       if (ititle.eq.0) write (iofile,91) >*/
  504.     if (ititle == 0) {
  505.     io__43.ciunit = status_1.iofile;
  506.     s_wsfe(&io__43);
  507.     e_wsfe();
  508.     }
  509. /*<       ititle=1 >*/
  510.     ititle = 1;
  511. /*<       write (iofile,afmt1) (anam(i),i=1,kntr) >*/
  512.     io__44.ciunit = status_1.iofile;
  513.     s_wsfe(&io__44);
  514.     i_1 = kntr;
  515.     for (i = 1; i <= i_1; ++i) {
  516.     do_fio(&c__1, (char *)&anam[i - 1], (ftnlen)sizeof(doublereal));
  517.     }
  518.     e_wsfe();
  519. /*<       write (iofile,afmt2) alstot,(vntot(i),i=1,kntr) >*/
  520.     io__46.ciunit = status_1.iofile;
  521.     s_wsfe(&io__46);
  522.     do_fio(&c__1, (char *)&alstot, (ftnlen)sizeof(doublereal));
  523.     i_1 = kntr;
  524.     for (i = 1; i <= i_1; ++i) {
  525.     do_fio(&c__1, (char *)&vntot[i - 1], (ftnlen)sizeof(doublereal));
  526.     }
  527.     e_wsfe();
  528. /*<   160 kntr=0 >*/
  529. L160:
  530.     kntr = 0;
  531. /*<       if (loc.ne.0) go to 120 >*/
  532.     if (loc != 0) {
  533.     goto L120;
  534.     }
  535.  
  536. /*  diodes */
  537.  
  538. /*<   200 if (jelcnt(11).eq.0) go to 300 >*/
  539. L200:
  540.     if (cirdat_1.jelcnt[10] == 0) {
  541.     goto L300;
  542.     }
  543. /*<       ititle=0 >*/
  544.     ititle = 0;
  545. /*<   201 format(//'0**** diode squared noise voltages (sq v/hz)') >*/
  546. /* L201: */
  547. /*<   210 loc=locate(11) >*/
  548. /* L210: */
  549.     loc = cirdat_1.locate[10];
  550. /*<       kntr=0 >*/
  551.     kntr = 0;
  552. /*<   220 if ((loc.eq.0).or.(nodplc(loc+16).ne.0)) go to 240 >*/
  553. L220:
  554.     if (loc == 0 || nodplc[loc + 15] != 0) {
  555.     goto L240;
  556.     }
  557. /*<       kntr=kntr+1 >*/
  558.     ++kntr;
  559. /*<       locv=nodplc(loc+1) >*/
  560.     locv = nodplc[loc];
  561. /*<       anam(kntr)=value(locv) >*/
  562.     anam[kntr - 1] = blank_1.value[locv - 1];
  563. /*<       node1=nodplc(loc+2) >*/
  564.     node1 = nodplc[loc + 1];
  565. /*<       node2=nodplc(loc+3) >*/
  566.     node2 = nodplc[loc + 2];
  567. /*<       node3=nodplc(loc+4) >*/
  568.     node3 = nodplc[loc + 3];
  569. /*<       locm=nodplc(loc+5) >*/
  570.     locm = nodplc[loc + 4];
  571. /*<       locm=nodplc(locm+1) >*/
  572.     locm = nodplc[locm];
  573. /*<       loct=nodplc(loc+11) >*/
  574.     loct = nodplc[loc + 10];
  575. /*<       area=value(locv+1) >*/
  576.     area = blank_1.value[locv];
  577. /*<       fnk=value(locm+10) >*/
  578.     fnk = blank_1.value[locm + 9];
  579. /*<       fna=value(locm+11) >*/
  580.     fna = blank_1.value[locm + 10];
  581.  
  582. /*  ohmic resistance */
  583.  
  584. /*<       cval=cvalue(lcvn+node1)-cvalue(lcvn+node3) >*/
  585.     i_1 = tabinf_1.lcvn + node1 - 1;
  586.     i_2 = tabinf_1.lcvn + node3 - 1;
  587.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  588.         .i;
  589.     cval->r = q_1.r, cval->i = q_1.i;
  590. /*<       vno1(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  591. /*<      1   *fourkt*value(locm+2)*area >*/
  592.     vno1[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  593.         *vimag)) * fourkt * blank_1.value[locm + 1] * area;
  594.  
  595. /*  junction shot noise and flicker noise */
  596.  
  597. /*<       cval=cvalue(lcvn+node3)-cvalue(lcvn+node2) >*/
  598.     i_1 = tabinf_1.lcvn + node3 - 1;
  599.     i_2 = tabinf_1.lcvn + node2 - 1;
  600.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  601.         .i;
  602.     cval->r = q_1.r, cval->i = q_1.i;
  603. /*<       vtemp=dble(vreal*vreal)+dble(vimag*vimag) >*/
  604.     vtemp = (doublereal) (*vreal * *vreal) + (doublereal) (*vimag * *vimag);
  605. /*<       arg=dmax1(dabs(value(lx0+loct+1)),1.0d-20) >*/
  606. /* Computing MAX */
  607.     d_2 = (d_1 = blank_1.value[tabinf_1.lx0 + loct], abs(d_1));
  608.     arg = max(1e-20,d_2);
  609. /*<       vno2(kntr)=vtemp*twoq*arg >*/
  610.     vno2[kntr - 1] = vtemp * twoq * arg;
  611. /*<       vno3(kntr)=vtemp*fnk*dexp(fna*dlog(arg))/freq >*/
  612.     vno3[kntr - 1] = vtemp * fnk * exp(fna * log(arg)) / freq;
  613. /*<       vntot(kntr)=vno1(kntr)+vno2(kntr)+vno3(kntr) >*/
  614.     vntot[kntr - 1] = vno1[kntr - 1] + vno2[kntr - 1] + vno3[kntr - 1];
  615. /*<       vnrms=vnrms+vntot(kntr) >*/
  616.     vnrms += vntot[kntr - 1];
  617. /*<       if (kntr.ge.kntlim) go to 250 >*/
  618.     if (kntr >= kntlim) {
  619.     goto L250;
  620.     }
  621. /*<   230 loc=nodplc(loc) >*/
  622. L230:
  623.     loc = nodplc[loc - 1];
  624. /*<       go to 220 >*/
  625.     goto L220;
  626. /*<   240 if (kntr.eq.0) go to 300 >*/
  627. L240:
  628.     if (kntr == 0) {
  629.     goto L300;
  630.     }
  631. /*<   250 if (nprnt.eq.0) go to 260 >*/
  632. L250:
  633.     if (nprnt == 0) {
  634.     goto L260;
  635.     }
  636. /*<       if (ititle.eq.0) write (iofile,201) >*/
  637.     if (ititle == 0) {
  638.     io__58.ciunit = status_1.iofile;
  639.     s_wsfe(&io__58);
  640.     e_wsfe();
  641.     }
  642. /*<       ititle=1 >*/
  643.     ititle = 1;
  644. /*<       write (iofile,afmt1) (anam(i),i=1,kntr) >*/
  645.     io__59.ciunit = status_1.iofile;
  646.     s_wsfe(&io__59);
  647.     i_1 = kntr;
  648.     for (i = 1; i <= i_1; ++i) {
  649.     do_fio(&c__1, (char *)&anam[i - 1], (ftnlen)sizeof(doublereal));
  650.     }
  651.     e_wsfe();
  652. /*<       write (iofile,afmt2) alsrs,(vno1(i),i=1,kntr) >*/
  653.     io__60.ciunit = status_1.iofile;
  654.     s_wsfe(&io__60);
  655.     do_fio(&c__1, (char *)&alsrs, (ftnlen)sizeof(doublereal));
  656.     i_1 = kntr;
  657.     for (i = 1; i <= i_1; ++i) {
  658.     do_fio(&c__1, (char *)&vno1[i - 1], (ftnlen)sizeof(doublereal));
  659.     }
  660.     e_wsfe();
  661. /*<       write (iofile,afmt2) alsid,(vno2(i),i=1,kntr) >*/
  662.     io__61.ciunit = status_1.iofile;
  663.     s_wsfe(&io__61);
  664.     do_fio(&c__1, (char *)&alsid, (ftnlen)sizeof(doublereal));
  665.     i_1 = kntr;
  666.     for (i = 1; i <= i_1; ++i) {
  667.     do_fio(&c__1, (char *)&vno2[i - 1], (ftnlen)sizeof(doublereal));
  668.     }
  669.     e_wsfe();
  670. /*<       write (iofile,afmt2) alsfn,(vno3(i),i=1,kntr) >*/
  671.     io__62.ciunit = status_1.iofile;
  672.     s_wsfe(&io__62);
  673.     do_fio(&c__1, (char *)&alsfn, (ftnlen)sizeof(doublereal));
  674.     i_1 = kntr;
  675.     for (i = 1; i <= i_1; ++i) {
  676.     do_fio(&c__1, (char *)&vno3[i - 1], (ftnlen)sizeof(doublereal));
  677.     }
  678.     e_wsfe();
  679. /*<       write (iofile,afmt2) alstot,(vntot(i),i=1,kntr) >*/
  680.     io__63.ciunit = status_1.iofile;
  681.     s_wsfe(&io__63);
  682.     do_fio(&c__1, (char *)&alstot, (ftnlen)sizeof(doublereal));
  683.     i_1 = kntr;
  684.     for (i = 1; i <= i_1; ++i) {
  685.     do_fio(&c__1, (char *)&vntot[i - 1], (ftnlen)sizeof(doublereal));
  686.     }
  687.     e_wsfe();
  688. /*<   260 kntr=0 >*/
  689. L260:
  690.     kntr = 0;
  691. /*<       if (loc.ne.0) go to 230 >*/
  692.     if (loc != 0) {
  693.     goto L230;
  694.     }
  695.  
  696. /*  bipolar junction transistors */
  697.  
  698. /*<   300 if (jelcnt(12).eq.0) go to 400 >*/
  699. L300:
  700.     if (cirdat_1.jelcnt[11] == 0) {
  701.     goto L400;
  702.     }
  703. /*<       ititle=0 >*/
  704.     ititle = 0;
  705. /*<   301 format(//'0**** transistor squared noise voltages (sq v/hz)') >*/
  706. /* L301: */
  707. /*<   310 loc=locate(12) >*/
  708. /* L310: */
  709.     loc = cirdat_1.locate[11];
  710. /*<       kntr=0 >*/
  711.     kntr = 0;
  712. /*<   320 if ((loc.eq.0).or.(nodplc(loc+36).ne.0)) go to 340 >*/
  713. L320:
  714.     if (loc == 0 || nodplc[loc + 35] != 0) {
  715.     goto L340;
  716.     }
  717. /*<       kntr=kntr+1 >*/
  718.     ++kntr;
  719. /*<       locv=nodplc(loc+1) >*/
  720.     locv = nodplc[loc];
  721. /*<       anam(kntr)=value(locv) >*/
  722.     anam[kntr - 1] = blank_1.value[locv - 1];
  723. /*<       node1=nodplc(loc+2) >*/
  724.     node1 = nodplc[loc + 1];
  725. /*<       node2=nodplc(loc+3) >*/
  726.     node2 = nodplc[loc + 2];
  727. /*<       node3=nodplc(loc+4) >*/
  728.     node3 = nodplc[loc + 3];
  729. /*<       node4=nodplc(loc+5) >*/
  730.     node4 = nodplc[loc + 4];
  731. /*<       node5=nodplc(loc+6) >*/
  732.     node5 = nodplc[loc + 5];
  733. /*<       node6=nodplc(loc+7) >*/
  734.     node6 = nodplc[loc + 6];
  735. /*<       locm=nodplc(loc+8) >*/
  736.     locm = nodplc[loc + 7];
  737. /*<       locm=nodplc(locm+1) >*/
  738.     locm = nodplc[locm];
  739. /*<       loct=nodplc(loc+22) >*/
  740.     loct = nodplc[loc + 21];
  741. /*<       area=value(locv+1) >*/
  742.     area = blank_1.value[locv];
  743. /*<       fnk=value(locm+44) >*/
  744.     fnk = blank_1.value[locm + 43];
  745. /*<       fna=value(locm+45) >*/
  746.     fna = blank_1.value[locm + 44];
  747.  
  748. /*  extrinsic resistances */
  749.  
  750. /* ...  base resistance */
  751. /*<       cval=cvalue(lcvn+node2)-cvalue(lcvn+node5) >*/
  752.     i_1 = tabinf_1.lcvn + node2 - 1;
  753.     i_2 = tabinf_1.lcvn + node5 - 1;
  754.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  755.         .i;
  756.     cval->r = q_1.r, cval->i = q_1.i;
  757. /*<       vno1(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  758. /*<      1  *fourkt*value(lx0+loct+16) >*/
  759.     vno1[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  760.         *vimag)) * fourkt * blank_1.value[tabinf_1.lx0 + loct + 15];
  761. /* ...  collector resistance */
  762. /*<       cval=cvalue(lcvn+node1)-cvalue(lcvn+node4) >*/
  763.     i_1 = tabinf_1.lcvn + node1 - 1;
  764.     i_2 = tabinf_1.lcvn + node4 - 1;
  765.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  766.         .i;
  767.     cval->r = q_1.r, cval->i = q_1.i;
  768. /*<       vno2(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  769. /*<      1   *fourkt*value(locm+20)*area >*/
  770.     vno2[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  771.         *vimag)) * fourkt * blank_1.value[locm + 19] * area;
  772. /* ...  emitter resistance */
  773. /*<       cval=cvalue(lcvn+node3)-cvalue(lcvn+node6) >*/
  774.     i_1 = tabinf_1.lcvn + node3 - 1;
  775.     i_2 = tabinf_1.lcvn + node6 - 1;
  776.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  777.         .i;
  778.     cval->r = q_1.r, cval->i = q_1.i;
  779. /*<       vno3(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  780. /*<      1   *fourkt*value(locm+19)*area >*/
  781.     vno3[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  782.         *vimag)) * fourkt * blank_1.value[locm + 18] * area;
  783.  
  784. /*  base current shot noise and flicker noise */
  785.  
  786. /*<       cval=cvalue(lcvn+node5)-cvalue(lcvn+node6) >*/
  787.     i_1 = tabinf_1.lcvn + node5 - 1;
  788.     i_2 = tabinf_1.lcvn + node6 - 1;
  789.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  790.         .i;
  791.     cval->r = q_1.r, cval->i = q_1.i;
  792. /*<       vtemp=dble(vreal*vreal)+dble(vimag*vimag) >*/
  793.     vtemp = (doublereal) (*vreal * *vreal) + (doublereal) (*vimag * *vimag);
  794. /*<       arg=dmax1(dabs(value(lx0+loct+3)),1.0d-20) >*/
  795. /* Computing MAX */
  796.     d_2 = (d_1 = blank_1.value[tabinf_1.lx0 + loct + 2], abs(d_1));
  797.     arg = max(1e-20,d_2);
  798. /*<       vno4(kntr)=vtemp*twoq*arg >*/
  799.     vno4[kntr - 1] = vtemp * twoq * arg;
  800. /*<       vno5(kntr)=vtemp*fnk*dexp(fna*dlog(arg))/freq >*/
  801.     vno5[kntr - 1] = vtemp * fnk * exp(fna * log(arg)) / freq;
  802.  
  803. /*  collector current shot noise */
  804.  
  805. /*<       cval=cvalue(lcvn+node4)-cvalue(lcvn+node6) >*/
  806.     i_1 = tabinf_1.lcvn + node4 - 1;
  807.     i_2 = tabinf_1.lcvn + node6 - 1;
  808.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  809.         .i;
  810.     cval->r = q_1.r, cval->i = q_1.i;
  811. /*<       vno6(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  812. /*<      1   *twoq*dabs(value(lx0+loct+2)) >*/
  813.     vno6[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  814.         *vimag)) * twoq * (d_1 = blank_1.value[tabinf_1.lx0 + loct + 1], 
  815.         abs(d_1));
  816. /*<       vntot(kntr)=vno1(kntr)+vno2(kntr)+vno3(kntr)+vno4(kntr)+vno5(kntr) >*/
  817. /*<      1   +vno6(kntr) >*/
  818.     vntot[kntr - 1] = vno1[kntr - 1] + vno2[kntr - 1] + vno3[kntr - 1] + vno4[
  819.         kntr - 1] + vno5[kntr - 1] + vno6[kntr - 1];
  820. /*<       vnrms=vnrms+vntot(kntr) >*/
  821.     vnrms += vntot[kntr - 1];
  822. /*<       if (kntr.ge.kntlim) go to 350 >*/
  823.     if (kntr >= kntlim) {
  824.     goto L350;
  825.     }
  826. /*<   330 loc=nodplc(loc) >*/
  827. L330:
  828.     loc = nodplc[loc - 1];
  829. /*<       go to 320 >*/
  830.     goto L320;
  831. /*<   340 if (kntr.eq.0) go to 400 >*/
  832. L340:
  833.     if (kntr == 0) {
  834.     goto L400;
  835.     }
  836. /*<   350 if (nprnt.eq.0) go to 360 >*/
  837. L350:
  838.     if (nprnt == 0) {
  839.     goto L360;
  840.     }
  841. /*<       if (ititle.eq.0) write (iofile,301) >*/
  842.     if (ititle == 0) {
  843.     io__70.ciunit = status_1.iofile;
  844.     s_wsfe(&io__70);
  845.     e_wsfe();
  846.     }
  847. /*<       ititle=1 >*/
  848.     ititle = 1;
  849. /*<       write (iofile,afmt1) (anam(i),i=1,kntr) >*/
  850.     io__71.ciunit = status_1.iofile;
  851.     s_wsfe(&io__71);
  852.     i_1 = kntr;
  853.     for (i = 1; i <= i_1; ++i) {
  854.     do_fio(&c__1, (char *)&anam[i - 1], (ftnlen)sizeof(doublereal));
  855.     }
  856.     e_wsfe();
  857. /*<       write (iofile,afmt2) alsrb,(vno1(i),i=1,kntr) >*/
  858.     io__72.ciunit = status_1.iofile;
  859.     s_wsfe(&io__72);
  860.     do_fio(&c__1, (char *)&alsrb, (ftnlen)sizeof(doublereal));
  861.     i_1 = kntr;
  862.     for (i = 1; i <= i_1; ++i) {
  863.     do_fio(&c__1, (char *)&vno1[i - 1], (ftnlen)sizeof(doublereal));
  864.     }
  865.     e_wsfe();
  866. /*<       write (iofile,afmt2) alsrc,(vno2(i),i=1,kntr) >*/
  867.     io__73.ciunit = status_1.iofile;
  868.     s_wsfe(&io__73);
  869.     do_fio(&c__1, (char *)&alsrc, (ftnlen)sizeof(doublereal));
  870.     i_1 = kntr;
  871.     for (i = 1; i <= i_1; ++i) {
  872.     do_fio(&c__1, (char *)&vno2[i - 1], (ftnlen)sizeof(doublereal));
  873.     }
  874.     e_wsfe();
  875. /*<       write (iofile,afmt2) alsre,(vno3(i),i=1,kntr) >*/
  876.     io__74.ciunit = status_1.iofile;
  877.     s_wsfe(&io__74);
  878.     do_fio(&c__1, (char *)&alsre, (ftnlen)sizeof(doublereal));
  879.     i_1 = kntr;
  880.     for (i = 1; i <= i_1; ++i) {
  881.     do_fio(&c__1, (char *)&vno3[i - 1], (ftnlen)sizeof(doublereal));
  882.     }
  883.     e_wsfe();
  884. /*<       write (iofile,afmt2) alsib,(vno4(i),i=1,kntr) >*/
  885.     io__75.ciunit = status_1.iofile;
  886.     s_wsfe(&io__75);
  887.     do_fio(&c__1, (char *)&alsib, (ftnlen)sizeof(doublereal));
  888.     i_1 = kntr;
  889.     for (i = 1; i <= i_1; ++i) {
  890.     do_fio(&c__1, (char *)&vno4[i - 1], (ftnlen)sizeof(doublereal));
  891.     }
  892.     e_wsfe();
  893. /*<       write (iofile,afmt2) alsic,(vno6(i),i=1,kntr) >*/
  894.     io__76.ciunit = status_1.iofile;
  895.     s_wsfe(&io__76);
  896.     do_fio(&c__1, (char *)&alsic, (ftnlen)sizeof(doublereal));
  897.     i_1 = kntr;
  898.     for (i = 1; i <= i_1; ++i) {
  899.     do_fio(&c__1, (char *)&vno6[i - 1], (ftnlen)sizeof(doublereal));
  900.     }
  901.     e_wsfe();
  902. /*<       write (iofile,afmt2) alsfn,(vno5(i),i=1,kntr) >*/
  903.     io__77.ciunit = status_1.iofile;
  904.     s_wsfe(&io__77);
  905.     do_fio(&c__1, (char *)&alsfn, (ftnlen)sizeof(doublereal));
  906.     i_1 = kntr;
  907.     for (i = 1; i <= i_1; ++i) {
  908.     do_fio(&c__1, (char *)&vno5[i - 1], (ftnlen)sizeof(doublereal));
  909.     }
  910.     e_wsfe();
  911. /*<       write (iofile,afmt2) alstot,(vntot(i),i=1,kntr) >*/
  912.     io__78.ciunit = status_1.iofile;
  913.     s_wsfe(&io__78);
  914.     do_fio(&c__1, (char *)&alstot, (ftnlen)sizeof(doublereal));
  915.     i_1 = kntr;
  916.     for (i = 1; i <= i_1; ++i) {
  917.     do_fio(&c__1, (char *)&vntot[i - 1], (ftnlen)sizeof(doublereal));
  918.     }
  919.     e_wsfe();
  920. /*<   360 kntr=0 >*/
  921. L360:
  922.     kntr = 0;
  923. /*<       if (loc.ne.0) go to 330 >*/
  924.     if (loc != 0) {
  925.     goto L330;
  926.     }
  927.  
  928. /*  jfets */
  929.  
  930. /*<   400 if (jelcnt(13).eq.0) go to 500 >*/
  931. L400:
  932.     if (cirdat_1.jelcnt[12] == 0) {
  933.     goto L500;
  934.     }
  935. /*<       ititle=0 >*/
  936.     ititle = 0;
  937. /*<   401 format(//'0**** jfet squared noise voltages (sq v/hz)') >*/
  938. /* L401: */
  939. /*<   410 loc=locate(13) >*/
  940. /* L410: */
  941.     loc = cirdat_1.locate[12];
  942. /*<       kntr=0 >*/
  943.     kntr = 0;
  944. /*<   420 if ((loc.eq.0).or.(nodplc(loc+25).ne.0)) go to 440 >*/
  945. L420:
  946.     if (loc == 0 || nodplc[loc + 24] != 0) {
  947.     goto L440;
  948.     }
  949. /*<       kntr=kntr+1 >*/
  950.     ++kntr;
  951. /*<       locv=nodplc(loc+1) >*/
  952.     locv = nodplc[loc];
  953. /*<       anam(kntr)=value(locv) >*/
  954.     anam[kntr - 1] = blank_1.value[locv - 1];
  955. /*<       node1=nodplc(loc+2) >*/
  956.     node1 = nodplc[loc + 1];
  957. /*<       node2=nodplc(loc+3) >*/
  958.     node2 = nodplc[loc + 2];
  959. /*<       node3=nodplc(loc+4) >*/
  960.     node3 = nodplc[loc + 3];
  961. /*<       node4=nodplc(loc+5) >*/
  962.     node4 = nodplc[loc + 4];
  963. /*<       node5=nodplc(loc+6) >*/
  964.     node5 = nodplc[loc + 5];
  965. /*<       locm=nodplc(loc+7) >*/
  966.     locm = nodplc[loc + 6];
  967. /*<       locm=nodplc(locm+1) >*/
  968.     locm = nodplc[locm];
  969. /*<       loct=nodplc(loc+19) >*/
  970.     loct = nodplc[loc + 18];
  971. /*<       area=value(locv+1) >*/
  972.     area = blank_1.value[locv];
  973. /*<       fnk=value(locm+10) >*/
  974.     fnk = blank_1.value[locm + 9];
  975. /*<       fna=value(locm+11) >*/
  976.     fna = blank_1.value[locm + 10];
  977.  
  978. /*  extrinsic resistances */
  979.  
  980. /* ...  drain resistance */
  981. /*<       cval=cvalue(lcvn+node1)-cvalue(lcvn+node4) >*/
  982.     i_1 = tabinf_1.lcvn + node1 - 1;
  983.     i_2 = tabinf_1.lcvn + node4 - 1;
  984.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  985.         .i;
  986.     cval->r = q_1.r, cval->i = q_1.i;
  987. /*<       vno1(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  988. /*<      1   *fourkt*value(locm+4)*area >*/
  989.     vno1[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  990.         *vimag)) * fourkt * blank_1.value[locm + 3] * area;
  991. /* ...  source resistance */
  992. /*<       cval=cvalue(lcvn+node3)-cvalue(lcvn+node5) >*/
  993.     i_1 = tabinf_1.lcvn + node3 - 1;
  994.     i_2 = tabinf_1.lcvn + node5 - 1;
  995.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  996.         .i;
  997.     cval->r = q_1.r, cval->i = q_1.i;
  998. /*<       vno2(kntr)=(dble(vreal*vreal)+dble(vimag*vimag)) >*/
  999. /*<      1   *fourkt*value(locm+5)*area >*/
  1000.     vno2[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  1001.         *vimag)) * fourkt * blank_1.value[locm + 4] * area;
  1002.  
  1003. /*  drain current shot noise and flicker noise */
  1004.  
  1005. /*<       cval=cvalue(lcvn+node4)-cvalue(lcvn+node5) >*/
  1006.     i_1 = tabinf_1.lcvn + node4 - 1;
  1007.     i_2 = tabinf_1.lcvn + node5 - 1;
  1008.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  1009.         .i;
  1010.     cval->r = q_1.r, cval->i = q_1.i;
  1011. /*<       vtemp=dble(vreal*vreal)+dble(vimag*vimag) >*/
  1012.     vtemp = (doublereal) (*vreal * *vreal) + (doublereal) (*vimag * *vimag);
  1013. /*<       vno3(kntr)=vtemp*fourkt*2.0d0*dabs(value(lx0+loct+5))/3.0d0 >*/
  1014.     vno3[kntr - 1] = vtemp * fourkt * 2. * (d_1 = blank_1.value[tabinf_1.lx0 
  1015.         + loct + 4], abs(d_1)) / 3.;
  1016. /*<       arg=dmax1(dabs(value(lx0+loct+3)),1.0d-20) >*/
  1017. /* Computing MAX */
  1018.     d_2 = (d_1 = blank_1.value[tabinf_1.lx0 + loct + 2], abs(d_1));
  1019.     arg = max(1e-20,d_2);
  1020. /*<       vno4(kntr)=vtemp*fnk*dexp(fna*dlog(arg))/freq >*/
  1021.     vno4[kntr - 1] = vtemp * fnk * exp(fna * log(arg)) / freq;
  1022. /*<       vntot(kntr)=vno1(kntr)+vno2(kntr)+vno3(kntr)+vno4(kntr) >*/
  1023.     vntot[kntr - 1] = vno1[kntr - 1] + vno2[kntr - 1] + vno3[kntr - 1] + vno4[
  1024.         kntr - 1];
  1025. /*<       vnrms=vnrms+vntot(kntr) >*/
  1026.     vnrms += vntot[kntr - 1];
  1027. /*<       if (kntr.ge.kntlim) go to 450 >*/
  1028.     if (kntr >= kntlim) {
  1029.     goto L450;
  1030.     }
  1031. /*<   430 loc=nodplc(loc) >*/
  1032. L430:
  1033.     loc = nodplc[loc - 1];
  1034. /*<       go to 420 >*/
  1035.     goto L420;
  1036. /*<   440 if (kntr.eq.0) go to 500 >*/
  1037. L440:
  1038.     if (kntr == 0) {
  1039.     goto L500;
  1040.     }
  1041. /*<   450 if (nprnt.eq.0) go to 460 >*/
  1042. L450:
  1043.     if (nprnt == 0) {
  1044.     goto L460;
  1045.     }
  1046. /*<       if (ititle.eq.0) write (iofile,401) >*/
  1047.     if (ititle == 0) {
  1048.     io__79.ciunit = status_1.iofile;
  1049.     s_wsfe(&io__79);
  1050.     e_wsfe();
  1051.     }
  1052. /*<       ititle=1 >*/
  1053.     ititle = 1;
  1054. /*<       write (iofile,afmt1) (anam(i),i=1,kntr) >*/
  1055.     io__80.ciunit = status_1.iofile;
  1056.     s_wsfe(&io__80);
  1057.     i_1 = kntr;
  1058.     for (i = 1; i <= i_1; ++i) {
  1059.     do_fio(&c__1, (char *)&anam[i - 1], (ftnlen)sizeof(doublereal));
  1060.     }
  1061.     e_wsfe();
  1062. /*<       write (iofile,afmt2) alsrd,(vno1(i),i=1,kntr) >*/
  1063.     io__81.ciunit = status_1.iofile;
  1064.     s_wsfe(&io__81);
  1065.     do_fio(&c__1, (char *)&alsrd, (ftnlen)sizeof(doublereal));
  1066.     i_1 = kntr;
  1067.     for (i = 1; i <= i_1; ++i) {
  1068.     do_fio(&c__1, (char *)&vno1[i - 1], (ftnlen)sizeof(doublereal));
  1069.     }
  1070.     e_wsfe();
  1071. /*<       write (iofile,afmt2) alsrs,(vno2(i),i=1,kntr) >*/
  1072.     io__82.ciunit = status_1.iofile;
  1073.     s_wsfe(&io__82);
  1074.     do_fio(&c__1, (char *)&alsrs, (ftnlen)sizeof(doublereal));
  1075.     i_1 = kntr;
  1076.     for (i = 1; i <= i_1; ++i) {
  1077.     do_fio(&c__1, (char *)&vno2[i - 1], (ftnlen)sizeof(doublereal));
  1078.     }
  1079.     e_wsfe();
  1080. /*<       write (iofile,afmt2) alsid,(vno3(i),i=1,kntr) >*/
  1081.     io__83.ciunit = status_1.iofile;
  1082.     s_wsfe(&io__83);
  1083.     do_fio(&c__1, (char *)&alsid, (ftnlen)sizeof(doublereal));
  1084.     i_1 = kntr;
  1085.     for (i = 1; i <= i_1; ++i) {
  1086.     do_fio(&c__1, (char *)&vno3[i - 1], (ftnlen)sizeof(doublereal));
  1087.     }
  1088.     e_wsfe();
  1089. /*<       write (iofile,afmt2) alsfn,(vno4(i),i=1,kntr) >*/
  1090.     io__84.ciunit = status_1.iofile;
  1091.     s_wsfe(&io__84);
  1092.     do_fio(&c__1, (char *)&alsfn, (ftnlen)sizeof(doublereal));
  1093.     i_1 = kntr;
  1094.     for (i = 1; i <= i_1; ++i) {
  1095.     do_fio(&c__1, (char *)&vno4[i - 1], (ftnlen)sizeof(doublereal));
  1096.     }
  1097.     e_wsfe();
  1098. /*<       write (iofile,afmt2) alstot,(vntot(i),i=1,kntr) >*/
  1099.     io__85.ciunit = status_1.iofile;
  1100.     s_wsfe(&io__85);
  1101.     do_fio(&c__1, (char *)&alstot, (ftnlen)sizeof(doublereal));
  1102.     i_1 = kntr;
  1103.     for (i = 1; i <= i_1; ++i) {
  1104.     do_fio(&c__1, (char *)&vntot[i - 1], (ftnlen)sizeof(doublereal));
  1105.     }
  1106.     e_wsfe();
  1107. /*<   460 kntr=0 >*/
  1108. L460:
  1109.     kntr = 0;
  1110. /*<       if (loc.ne.0) go to 430 >*/
  1111.     if (loc != 0) {
  1112.     goto L430;
  1113.     }
  1114.  
  1115. /*  mosfets */
  1116.  
  1117. /*<   500 if (jelcnt(14).eq.0) go to 600 >*/
  1118. L500:
  1119.     if (cirdat_1.jelcnt[13] == 0) {
  1120.     goto L600;
  1121.     }
  1122. /*<       ititle=0 >*/
  1123.     ititle = 0;
  1124. /*<   501 format(//'0**** mosfet squared noise voltages (sq v/hz)') >*/
  1125. /* L501: */
  1126. /*<   510 loc=locate(14) >*/
  1127. /* L510: */
  1128.     loc = cirdat_1.locate[13];
  1129. /*<       kntr=0 >*/
  1130.     kntr = 0;
  1131. /*<   520 if ((loc.eq.0).or.(nodplc(loc+33).ne.0)) go to 540 >*/
  1132. L520:
  1133.     if (loc == 0 || nodplc[loc + 32] != 0) {
  1134.     goto L540;
  1135.     }
  1136. /*<       kntr=kntr+1 >*/
  1137.     ++kntr;
  1138. /*<       locv=nodplc(loc+1) >*/
  1139.     locv = nodplc[loc];
  1140. /*<       anam(kntr)=value(locv) >*/
  1141.     anam[kntr - 1] = blank_1.value[locv - 1];
  1142. /*<       node1=nodplc(loc+2) >*/
  1143.     node1 = nodplc[loc + 1];
  1144. /*<       node2=nodplc(loc+3) >*/
  1145.     node2 = nodplc[loc + 2];
  1146. /*<       node3=nodplc(loc+4) >*/
  1147.     node3 = nodplc[loc + 3];
  1148. /*<       node4=nodplc(loc+5) >*/
  1149.     node4 = nodplc[loc + 4];
  1150. /*<       node5=nodplc(loc+6) >*/
  1151.     node5 = nodplc[loc + 5];
  1152. /*<       node6=nodplc(loc+7) >*/
  1153.     node6 = nodplc[loc + 6];
  1154. /*<       locm=nodplc(loc+8) >*/
  1155.     locm = nodplc[loc + 7];
  1156. /*<       itype=nodplc(locm+2) >*/
  1157.     itype = nodplc[locm + 1];
  1158. /*<       loct=nodplc(loc+26) >*/
  1159.     loct = nodplc[loc + 25];
  1160. /*<       locm=nodplc(locm+1) >*/
  1161.     locm = nodplc[locm];
  1162. /*<       xl=value(locv+1)-2.0d0*value(locm+28) >*/
  1163.     xl = blank_1.value[locv] - blank_1.value[locm + 27] * 2.;
  1164. /*<       xw=value(locv+2) >*/
  1165.     xw = blank_1.value[locv + 1];
  1166. /*<       cox=value(locm+22) >*/
  1167.     cox = blank_1.value[locm + 21];
  1168. /*<       if (cox.le.0.0d0) cox=epsox/1.0d-7 >*/
  1169.     if (cox <= 0.) {
  1170.     cox = knstnt_1.epsox / 1e-7;
  1171.     }
  1172. /*<       fnk=value(locm+36) >*/
  1173.     fnk = blank_1.value[locm + 35];
  1174. /*<       fna=value(locm+37) >*/
  1175.     fna = blank_1.value[locm + 36];
  1176.  
  1177. /*  extrinsic resistances */
  1178.  
  1179. /*<       if ((value(locm+7).le.0.0d0).and. >*/
  1180. /*<      1   (value(locm+8).le.0.0d0)) go to 522 >*/
  1181.     if (blank_1.value[locm + 6] <= 0. && blank_1.value[locm + 7] <= 0.) {
  1182.     goto L522;
  1183.     }
  1184. /*<       gdpr=value(locm+7) >*/
  1185.     gdpr = blank_1.value[locm + 6];
  1186. /*<       gspr=value(locm+8) >*/
  1187.     gspr = blank_1.value[locm + 7];
  1188. /*<       go to 524 >*/
  1189.     goto L524;
  1190. /*<   522 gdpr=value(locm+16)/value(locv+13) >*/
  1191. L522:
  1192.     gdpr = blank_1.value[locm + 15] / blank_1.value[locv + 12];
  1193. /*<       gspr=value(locm+16)/value(locv+14) >*/
  1194.     gspr = blank_1.value[locm + 15] / blank_1.value[locv + 13];
  1195. /* ...  drain resistance */
  1196. /*<   524 cval=cvalue(lcvn+node1)-cvalue(lcvn+node5) >*/
  1197. L524:
  1198.     i_1 = tabinf_1.lcvn + node1 - 1;
  1199.     i_2 = tabinf_1.lcvn + node5 - 1;
  1200.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  1201.         .i;
  1202.     cval->r = q_1.r, cval->i = q_1.i;
  1203. /*<       vno1(kntr)=(dble(vreal*vreal)+dble(vimag*vimag))*fourkt*gdpr >*/
  1204.     vno1[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  1205.         *vimag)) * fourkt * gdpr;
  1206. /* ...  source resistance */
  1207. /*<       cval=cvalue(lcvn+node3)-cvalue(lcvn+node6) >*/
  1208.     i_1 = tabinf_1.lcvn + node3 - 1;
  1209.     i_2 = tabinf_1.lcvn + node6 - 1;
  1210.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  1211.         .i;
  1212.     cval->r = q_1.r, cval->i = q_1.i;
  1213. /*<       vno2(kntr)=(dble(vreal*vreal)+dble(vimag*vimag))*fourkt*gspr >*/
  1214.     vno2[kntr - 1] = ((doublereal) (*vreal * *vreal) + (doublereal) (*vimag * 
  1215.         *vimag)) * fourkt * gspr;
  1216.  
  1217. /*  drain current shot noise and flicker noise */
  1218.  
  1219. /*<       cval=cvalue(lcvn+node5)-cvalue(lcvn+node6) >*/
  1220.     i_1 = tabinf_1.lcvn + node5 - 1;
  1221.     i_2 = tabinf_1.lcvn + node6 - 1;
  1222.     q_1.r = cvalue[i_1].r - cvalue[i_2].r, q_1.i = cvalue[i_1].i - cvalue[i_2]
  1223.         .i;
  1224.     cval->r = q_1.r, cval->i = q_1.i;
  1225. /*<       vtemp=dble(vreal*vreal)+dble(vimag*vimag) >*/
  1226.     vtemp = (doublereal) (*vreal * *vreal) + (doublereal) (*vimag * *vimag);
  1227. /*<       gm=value(lx0+loct+7) >*/
  1228.     gm = blank_1.value[tabinf_1.lx0 + loct + 6];
  1229. /*<       arg=dmax1(dabs(value(lx0+loct+4)),1.0d-20) >*/
  1230. /* Computing MAX */
  1231.     d_2 = (d_1 = blank_1.value[tabinf_1.lx0 + loct + 3], abs(d_1));
  1232.     arg = max(1e-20,d_2);
  1233. /*<       vno3(kntr)=vtemp*fourkt*dabs(gm)/1.5d0 >*/
  1234.     vno3[kntr - 1] = vtemp * fourkt * abs(gm) / 1.5;
  1235. /*<       vno4(kntr)=vtemp*fnk*dexp(fna*dlog(arg))/(freq*cox*xl*xl) >*/
  1236.     vno4[kntr - 1] = vtemp * fnk * exp(fna * log(arg)) / (freq * cox * xl * 
  1237.         xl);
  1238. /*<   525 vntot(kntr)=vno1(kntr)+vno2(kntr)+vno3(kntr)+vno4(kntr) >*/
  1239. /* L525: */
  1240.     vntot[kntr - 1] = vno1[kntr - 1] + vno2[kntr - 1] + vno3[kntr - 1] + vno4[
  1241.         kntr - 1];
  1242. /*<       vnrms=vnrms+vntot(kntr) >*/
  1243.     vnrms += vntot[kntr - 1];
  1244. /*<       if (kntr.ge.kntlim) go to 550 >*/
  1245.     if (kntr >= kntlim) {
  1246.     goto L550;
  1247.     }
  1248. /*<   530 loc=nodplc(loc) >*/
  1249. L530:
  1250.     loc = nodplc[loc - 1];
  1251. /*<       go to 520 >*/
  1252.     goto L520;
  1253. /*<   540 if (kntr.eq.0) go to 600 >*/
  1254. L540:
  1255.     if (kntr == 0) {
  1256.     goto L600;
  1257.     }
  1258. /*<   550 if (nprnt.eq.0) go to 560 >*/
  1259. L550:
  1260.     if (nprnt == 0) {
  1261.     goto L560;
  1262.     }
  1263. /*<       if (ititle.eq.0) write (iofile,501) >*/
  1264.     if (ititle == 0) {
  1265.     io__93.ciunit = status_1.iofile;
  1266.     s_wsfe(&io__93);
  1267.     e_wsfe();
  1268.     }
  1269. /*<       ititle=1 >*/
  1270.     ititle = 1;
  1271. /*<       write (iofile,afmt1) (anam(i),i=1,kntr) >*/
  1272.     io__94.ciunit = status_1.iofile;
  1273.     s_wsfe(&io__94);
  1274.     i_1 = kntr;
  1275.     for (i = 1; i <= i_1; ++i) {
  1276.     do_fio(&c__1, (char *)&anam[i - 1], (ftnlen)sizeof(doublereal));
  1277.     }
  1278.     e_wsfe();
  1279. /*<       write (iofile,afmt2) alsrd,(vno1(i),i=1,kntr) >*/
  1280.     io__95.ciunit = status_1.iofile;
  1281.     s_wsfe(&io__95);
  1282.     do_fio(&c__1, (char *)&alsrd, (ftnlen)sizeof(doublereal));
  1283.     i_1 = kntr;
  1284.     for (i = 1; i <= i_1; ++i) {
  1285.     do_fio(&c__1, (char *)&vno1[i - 1], (ftnlen)sizeof(doublereal));
  1286.     }
  1287.     e_wsfe();
  1288. /*<       write (iofile,afmt2) alsrs,(vno2(i),i=1,kntr) >*/
  1289.     io__96.ciunit = status_1.iofile;
  1290.     s_wsfe(&io__96);
  1291.     do_fio(&c__1, (char *)&alsrs, (ftnlen)sizeof(doublereal));
  1292.     i_1 = kntr;
  1293.     for (i = 1; i <= i_1; ++i) {
  1294.     do_fio(&c__1, (char *)&vno2[i - 1], (ftnlen)sizeof(doublereal));
  1295.     }
  1296.     e_wsfe();
  1297. /*<       write (iofile,afmt2) alsid,(vno3(i),i=1,kntr) >*/
  1298.     io__97.ciunit = status_1.iofile;
  1299.     s_wsfe(&io__97);
  1300.     do_fio(&c__1, (char *)&alsid, (ftnlen)sizeof(doublereal));
  1301.     i_1 = kntr;
  1302.     for (i = 1; i <= i_1; ++i) {
  1303.     do_fio(&c__1, (char *)&vno3[i - 1], (ftnlen)sizeof(doublereal));
  1304.     }
  1305.     e_wsfe();
  1306. /*<       write (iofile,afmt2) alsfn,(vno4(i),i=1,kntr) >*/
  1307.     io__98.ciunit = status_1.iofile;
  1308.     s_wsfe(&io__98);
  1309.     do_fio(&c__1, (char *)&alsfn, (ftnlen)sizeof(doublereal));
  1310.     i_1 = kntr;
  1311.     for (i = 1; i <= i_1; ++i) {
  1312.     do_fio(&c__1, (char *)&vno4[i - 1], (ftnlen)sizeof(doublereal));
  1313.     }
  1314.     e_wsfe();
  1315. /*<       write (iofile,afmt2) alstot,(vntot(i),i=1,kntr) >*/
  1316.     io__99.ciunit = status_1.iofile;
  1317.     s_wsfe(&io__99);
  1318.     do_fio(&c__1, (char *)&alstot, (ftnlen)sizeof(doublereal));
  1319.     i_1 = kntr;
  1320.     for (i = 1; i <= i_1; ++i) {
  1321.     do_fio(&c__1, (char *)&vntot[i - 1], (ftnlen)sizeof(doublereal));
  1322.     }
  1323.     e_wsfe();
  1324. /*<   560 kntr=0 >*/
  1325. L560:
  1326.     kntr = 0;
  1327. /*<       if (loc.ne.0) go to 530 >*/
  1328.     if (loc != 0) {
  1329.     goto L530;
  1330.     }
  1331.  
  1332. /*  compute equivalent input noise voltage */
  1333.  
  1334. /*<   600 vnout=dsqrt(vnrms) >*/
  1335. L600:
  1336.     vnout = sqrt(vnrms);
  1337. /*<       vnin=vnout/vout >*/
  1338.     vnin = vnout / vout;
  1339. /*<       if (nprnt.eq.0) go to 620 >*/
  1340.     if (nprnt == 0) {
  1341.     goto L620;
  1342.     }
  1343. /*<       do 610 i=1,5 >*/
  1344.     for (i = 1; i <= 5; ++i) {
  1345. /*<       string(i)=ablnk >*/
  1346.     string[i - 1] = ablnk;
  1347. /*<   610 continue >*/
  1348. /* L610: */
  1349.     }
  1350. /*<       ioutyp=1 >*/
  1351.     ioutyp = 1;
  1352. /*<       ipos=1 >*/
  1353.     ipos = 1;
  1354. /*<       call outnam(nosout,ioutyp,string,ipos) >*/
  1355.     outnam_(&ac_1.nosout, &ioutyp, string, &ipos);
  1356. /*<       call move(string,ipos,aslash,1,1) >*/
  1357.     move_(string, &ipos, &aslash, &c__1, &c__1);
  1358. /*<       ipos=ipos+1 >*/
  1359.     ++ipos;
  1360. /*<       locv=nodplc(nosin+1) >*/
  1361.     locv = nodplc[ac_1.nosin];
  1362. /*<       anam1=value(locv) >*/
  1363.     anam1 = blank_1.value[locv - 1];
  1364. /*<       call move(string,ipos,anam1,1,8) >*/
  1365.     move_(string, &ipos, &anam1, &c__1, &c__8);
  1366. /*<       write (iofile,611) vnrms,vnout,string,vout,anam1,vnin >*/
  1367.     io__105.ciunit = status_1.iofile;
  1368.     s_wsfe(&io__105);
  1369.     do_fio(&c__1, (char *)&vnrms, (ftnlen)sizeof(doublereal));
  1370.     do_fio(&c__1, (char *)&vnout, (ftnlen)sizeof(doublereal));
  1371.     do_fio(&c__5, (char *)&string[0], (ftnlen)sizeof(doublereal));
  1372.     do_fio(&c__1, (char *)&vout, (ftnlen)sizeof(doublereal));
  1373.     do_fio(&c__1, (char *)&anam1, (ftnlen)sizeof(doublereal));
  1374.     do_fio(&c__1, (char *)&vnin, (ftnlen)sizeof(doublereal));
  1375.     e_wsfe();
  1376. /*<   611 format(////, >*/
  1377. /*<      1   '0**** total output noise voltage',9x,'= ',1pd10.3,' sq v/hz'/, >*/
  1378. /*<      2   1h0,40x,'= ',d10.3,' v/rt hz'/, >*/
  1379. /*<      3   '0     transfer function value:',/, >*/
  1380. /*<      4   1h0,7x,4a8,a1,'= ',d10.3,/, >*/
  1381. /*<      5   '0     equivalent input noise at ',a8,' = ',d10.3,' /rt hz') >*/
  1382.  
  1383. /*  save noise outputs */
  1384.  
  1385. /*<   620 loc=locate(44) >*/
  1386. L620:
  1387.     loc = cirdat_1.locate[43];
  1388. /*<   630 if (loc.eq.0) go to 1000 >*/
  1389. L630:
  1390.     if (loc == 0) {
  1391.     goto L1000;
  1392.     }
  1393. /*<       iseq=nodplc(loc+4) >*/
  1394.     tabinf_1.iseq = nodplc[loc + 3];
  1395. /*<       if (nodplc(loc+5).ne.2) go to 640 >*/
  1396.     if (nodplc[loc + 4] != 2) {
  1397.     goto L640;
  1398.     }
  1399. /*<       cvalue(loco+iseq)=vnout >*/
  1400.     i_1 = *loco + tabinf_1.iseq - 1;
  1401.     cvalue[i_1].r = vnout, cvalue[i_1].i = (float)0.;
  1402. /*<       go to 650 >*/
  1403.     goto L650;
  1404. /*<   640 cvalue(loco+iseq)=vnin >*/
  1405. L640:
  1406.     i_1 = *loco + tabinf_1.iseq - 1;
  1407.     cvalue[i_1].r = vnin, cvalue[i_1].i = (float)0.;
  1408. /*<   650 loc=nodplc(loc) >*/
  1409. L650:
  1410.     loc = nodplc[loc - 1];
  1411. /*<       go to 630 >*/
  1412.     goto L630;
  1413.  
  1414. /*  finished */
  1415.  
  1416. /*<  1000 return >*/
  1417. L1000:
  1418.     return 0;
  1419. /*<       end >*/
  1420. } /* noise_ */
  1421.  
  1422. #undef cvalue
  1423. #undef nodplc
  1424. #undef vreal
  1425. #undef vimag
  1426. #undef v
  1427. #undef c
  1428. #undef cval
  1429. #undef alsrb
  1430. #undef titln
  1431. #undef afmt2
  1432. #undef afmt1
  1433. #undef ablnk
  1434. #undef aslash
  1435. #undef alstot
  1436. #undef alsfn
  1437. #undef alsid
  1438. #undef alsic
  1439. #undef alsib
  1440. #undef alsrd
  1441. #undef alsrs
  1442. #undef alsre
  1443. #undef alsrc
  1444.  
  1445.  
  1446.